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ABSTRACT 


An  atom  adsorbed  on  the  surface  of  a  harmonic  crystal  is  considered.  The 


binding  potential  well  supports  many  vibrational  bound-states,  and  the  coupling 


of  the  motion  of  the  atom  with  the  substrate  gives  rise  to  phonon -exchange 


reactions,  which  subsequently  amount  to  atomic  transitions  between  the  bound 


states.  This  process  of  thermal  relaxation  of  the  adatom  density  operator  is 


commonly  described  with  reservoir  theory,  in  which  the  crystal  is  regarded  as  a 


thermal  bath  with  an  extremely  short  correlation  time.  The  latter  property  then 


justifies  the  neglect  of  any  memory  in  the  interaction,  which  is  usually  referred 


to  as  the  Markov  approximation.  It  is  shown,  however,  that  the  reservoir 


correlation  time  is  not  necessarily  small  in  comparison  with  the  inverse 


(relaxation)  frequencies  of  the  system,  which  implies  a  breakdown  of  this 


approach.  The  equation  of  motion  for  the  adatom  density-operator  is  then  solved 


without  the  Markov  approximation,  and  the  result  is  used  for  the  evaluation  of 


the  spectrum  and  line  profile  for  absorption  of  infrared  radiation.  Essentially 


different  results  are  obtained  than  with  the  approximate  theory,  but  it  turns  out 


that  the  solutions  can  be  expressed  in  terms  of  the  same  parameters. 
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I.  INTRODUCTION 


A  very  sensitive  method  for  studying  atomic  or  molecular  systems  and  the 
interaction  with  their  environment  (for  instance  collisions)  is  by  observation  of 
the  absorption  line  shapes.  A  low-intensity  laser  with  power  1^  (energy  per  unit 
of  time  through  a  unit  area,  perpendicular  to  the  direction  of  propagation). 


frequency  and  polarization  is  scanned  over  the  resonances  of  the  molecule. 
Dipole  coupling  between  the  system  and  the  laser  field  then  amounts  to  absorption 
of  radiation  at  a  rate  (energy  per  unit  of  time)^ 


I(uL) 


•  iu ,T  * 

dt  e  Trp{w(T)-eL,u-eLl 


(1.1) 


with  w(t)  the  dipole -moment  operator  in  the  Heisenberg  picture.  The  thermal- 
equilibrium  density  operator  p  represents  the  state  of  the  entire  system  of  the 
molecule,  environment  and  interaction,  but  not  the  laser.  The  advantage  of  weak- 
field  absorption  is  that  the  incident  field  does  not  disturb  the  system,  but  only 
probes  it.  Then  the  absorption  profile  can  be  obtained  from  the  Golden  Rule 
which  leads  to  Eq.  (1.1),  where  reference  to  the  laser  field  in  the  dynamical 
variables  p  and  y(x)  has  disappeared. 

Ue  consider  adsorbed  atoms  on  the  surface  of  a  crystal.  Electromagnetic 
(van  der  Uaals)  interaction  between  an  adsorbate  and  the  atoms  of  the  substrate 
is  accounted  for  by  a  potential  well,  which  supports  many  bounds  states.  The 
motion  of  the  atom  will  be  restricted  mainly  to  the  direction  perpendicular  to 
the  surface,  which  will  be  denoted  by  the  unit  vector  e  .  Transition  frequencies 
between  vibrational  states  are  in  the  infrared,  and  the  motion  of  the  atom 
induces  an  optical  activity  in  this  region  of  the  spectrum.  This  implies  that 
the  bond  has  non-vanishing  dipole-moment  matrix  elements  between  the  various 
states,  whereas  both  the  atom  and  the  crystal  are  assumed  to  be  transparent  for 
infrared  light.  Since  we  neglect  lateral  motion,  the  dipole-moment  operator 


takes  the  fora  y  ■  ye  ,  with  y  a  scalar  operator  in  the  Hilbert  space,  which  is 

*■*  ^  z 

spanned  by  the  adatom  states.  Then  we  can  write 

y(x)  -  e  exp(iLx)y  ,  (1.2) 

""  ""  Z 


for  the  time  evolution  of  y,  where  L  is  the  Liouvillian,  which  is  related  to  the 
Hamiltonian  according  to 

Lp  -  *_1[H,p]  ,  (1.3) 


for  an  arbitrary  Hilbert-space  operator  p.  Here,  H  represents  the  entire  system 
of  atom,  crystal  and  interaction.  Subsequently  we  use  the  identity 

Trp(y(t),y]  ■  Tryexp(-iLx){y,p]  ,  (1.4) 

which  allows  to  cast  the  absorption  profile  in  the  fora 


[(uL) 


*  I.H 


IV*J‘ 


L  L  e  he 
o 


p 

Jn 


Kid, -L)t 
Re  Tr  y  I  dx  e  (y.p) 

'0 


(1.5) 


This  representation  clearly  shows  that  the  peaks  of  the  spectrum  are  located  at 
the  resonances  of  the  complete  system,  which  equal  the  eigenvalues  of  L  rather 
than  the  adatom  transition  frequencies.  Moreover,  the  frequency  dependence  of 
ICu^)  reveals  the  details  of  the  dynamics.  From  Eq.  (1.5)  it  follows  that  the 
shape  of  a  spectral  line  is  germane  to  the  time  evolution  of  the  density 
operator,  which  explains  the  significance  of  the  study  of  absorption  spectra. 


II.  EQUATION  OF  MOTION 

The  atom  with  mass  m  is  bounded  to  the  crystal  by  a  potential  V(z),  where  z 
is  the  atom-surface  normal  distance,  and  hence  the  adbond  Hamiltonian  takes  the 
explicit  fora 


We  shall  not  specify  V(z)  any  further,  but  only  remark  that  the  common  choice  is 
2-7 

a  Norse  potential.  A  general  representation  for  the  Hamiltonian  of  a  harmonic 
crystal  reads** 


S-I 


(2.2) 


in  terms  of  the  phonon  creation  (a£g)  and  annihilation  (a^g)  operators,  which 
obey  boson  commutation  relations.  The  summation  runs  over  the  modes  ks,  with  k 
the  wave  vector  and  s  the  polarization,  and  wg(k)  is  the  dispersion  relation. 
Coupling  between  the  adatom  bond  and  the  phonon  field  is  assumed  to  be  brought 
about  by  single-phonon  transitions,  which  is  accurate,  as  long  as  the  level 
separations  of  Hg  do  not  exceed  the  cut-off  frequency  of  the  dispersion  relation. 
Then  the  complete  Hamiltonian  can  be  written  as 


H  -  H  +  H  -  RS 

«  P 


(2.3) 


with  the  atomic  part  of  the  coupling  given  by' 


S  -  dV/dz 


(2.4) 


in  terms  of  the  prescribed  potential  V(z).  The  operator  R  equals  the  z-component 
of  the  displacement  field  for  the  crystal  atoms  evaluated  in  the  vicinity  of  the 
adsorbate.  Explicitly  this  becomes 


2MVug(k)  (*ks  +  *-ks)(?ks'?t)  * 


(2.5) 


which  involves  the  mass  M  of  a  crystal  atom,  the  volumes  V  and  v  of  the  crystal 


V  la  m."  a.*  a_"  aJ*  a/  mJ? a  *  a*  l'.  .  i  .  »  . 


and  a  unit  cell,  respectively,  and  the  unit  polarization  vectors  e^. 

The  equation  of  motion  for  the  density  operator  p(t)  is 

il^p(t)  -  [H,p(t)]  ,  (2.6) 

which  has  the  solution 

p(t)  ■  exp(-iLt)p(0)  (2.7) 

in  terms  of  a  given  intial  state  p(0).  Comparison  with  expression  (1.5)  then 
shows  that  the  time  evolution  of  p(t)  is  governed  by  the  same  exponential  which 
determines  the  spectral  distribution  of  the  absorption.  Futhermore,  Eq.  (1.5) 
contains  the  thermal-equilibrium  state  p,  which  equals  the  long-time  solution 

p  *  lim  p(t)  .  (2.8) 

t-**> 

HI.  RESERVOIR  INTEGRAL 

It  will  be  obvious  that  an  exact  evaluation  of  exp(-iLt)  is  not  feasible. 
Fortunately,  this  is  not  necessary.  It  will  turn  out  to  be  sufficient  for  the 
calcualtion  of  observable  quantities,  like  I(u^),  to  obtain  an  equation  for  the 
reduced  adatom  density-operator,  defined  by 

P.(t)  ■  Tr  p(t)  ,  (3.1) 

o  p 

where  the  trace  runs  over  the  states  of  the  phonon  field.  The  large  crystal  with 
its  many,  closely-spaced,  degrees  of  freedom  (phonon  modes)  merely  acts  as  a 
thermal  bath,  and  its  thermal-equilibrium  density  operator  p^  is  not  affected  by 
the  presence  of  the  single  atom. 

9-11 

Now  it  is  a  standard  procedure  in  reservoir  theory  to  derive  an  integral 
of  the  equation  of  motion  (2.6),  which  only  involves  the  adbond  density  operator 
Po(t).  In  particular,  the  product-form  -RS  for  the  interaction  allows  a  concise 


6 


formulation,  and  the  result  is1*  (see  Appendix  for  details) 

i  4r  p  (t)  ■  L  p  (t)  ■  L_  [dr  exp(-iL  x)L  (x)p  (t-x)  ,  <  3 . i 

at  o  ao  tit  S  ac  o 


where  we  have  introduced  the  Liouvillians 


Lfip  -  fi  [Ha,p] 


LgP  ■  [S,p] 


(3.3) 


(3.4) 


pertaining  to  the  free  evolution  of  the  adatom  (Lfl)  and  the  atomic  part  of  the 

interaction  (S).  The  Liouvillian  L  (x)  is  defined  as 

c 


L  (x)p  «  G(x)Sp  -  G(x)  pS  , 
c 


which  contains  the  reservoir  correlation  function 


(3.5) 


G(x)  »  2irh  Tr  Rp  exp(iL  x)R  , 
P  P  P 


(3.6) 


with  LpP  *  H  [Hp,p].  Properties  of  the  substrate  only  enter  the  equation  (3.2) 
through  the  function  G(x),  which  equals  the  correlation  function  of  the 
displacement  of  the  surface,  due  to  thermal  motion  of  the  crystal  atoms.  In 
Section  VIII  we  shall  evaluate  G(x)  explicitly  for  a  simple  model. 

In  deriving  Eq.  (3.2)  from  Eq.  (2.6)  we  have  taken  for  the  initial  state  the 
factorized  form 


P(0)  ■  p  (0)5  , 

o  p 


(3.7) 


which  can  be  done  arbitrarily,  since  p(0)  is  not  prescribed  by  the  equation  of 
motion.  The  only  approximation  which  had  to  be  made  to  arrive  at  Eq.  (3.2)  was  a 
factorization  p(t-x)  -  pQ(t-x)Pp  in  the  integrand.  Of  course,  this  is  not  exact, 
but  it  is  much  better  than  the  usual  Markov  approximation  (Section  VI),  which  can 


'*»»» 


’  -  ■/  *  •  V  .•  V  V  V  V  V  'j-  ’ 


be  argued  as  follows.  The  first  term  on  the  right-hand  side  of  Eq.  (3.2), 

L^o^t),  would  account  for  the  evolution  of  the  adatom  if  there  were  not  crystal 

at  all  (RS  =  0),  and  the  second  term  represents  the  coupling  to  the  heat  bath. 

This  interaction  gives  rise  to  phonon  transition  between  levels  |k>  and  |t>, 

which  can  be  regarded  as  stimulated  absorptions  of  phonons  from  the  crystal 

(excitation  of  the  atomic  bond),  and  emissions  of  phonons  into  the  crystal  (decay 

of  the  adsorbate).  A  transition  Jk>  -*■  |l>  then  occurs  at  a  rate  n^a^,  with  nfc 

the  population  of  level  |k>  and  a^  the  rate  constant  for  this  particular 

transition.  If  we  look  at  the  relaxation  integral  in  Eq.  (3.2)  as  an  operator 

acting  on  pQ(t-T),  then  a^  equals  a  matrix  element  of  this  operator.  The 

factorization  p(t-x)  -  p  (t-x)p  in  the  integrand  then  introduces  a  small  error 

-  o  p 

in  a^*  but  the  correct  time  dependence  is  retained.  Because  it  is  the  time 
evolution  of  the  density  operator  which  determines  the  spectral  distribution  of 
the  absorption,  rather  than  the  density  operator  itself,  only  a  small  err^r  is 
made  by  the  factorization. 


IV.  SOLUTION 

In  this  section  we  solve  Eq.  (3.2).  To  this  end,  we  first  introduce  the 
Fourier-Laplace  transform  of  G(t)  by 

G(w)  ■  x  *  [  dt  exp(iwx)G(x)  ,  (4.1) 

J0 

which  has  the  inverse  integral 


and  in  the  same  way  we  define  the  Fourier-Laplace  transform  of  any  other  time- 
dependent  quantity.  Substitution  of  the  inverse  integral  of  Lc(t)  into  the 
equation  of  motion  (3.2)  then  yields 


in  terms  of  an  operator  inversion.  Here,  the  frequency-dependent  relaxation 
operator  f(u)  is  defined  as 


SLsf 
J  -01 


i 

oj-oj'  -L 


L  (w1  ) 

c 


(4.6) 


Equation  (4.5)  relates  the  (transformed)  adatom  density  operator  for  t  >  0  to  the 

given  initial  value  p  (0),  which  is  the  desired  solution. 

o 

Although  Eq.  (4.5)  is  an  explicit  solution,  the  combination  of  the  operator 
inversion,  with  the  definition  of  r(w)  as  an  integral  over  a  product  of  an 
inverse  operator  with  Lc(u),  might  seem  awkward.  We  shall  show,  however,  that 
the  evaluation  of  the  absorption  profile  has  exactly  the  same  degree  of 
complexity  as  in  the  usual  Markov  approximation. 


V.  EVALUATION  OF  r(u>) 

A  matrix  representation  of  T(u)  is  easily  derived.  First  we  rewite  Eq. 
(4.6)  in  the  form 


r<“)  ■  h  h  jV  r  * 


i((i)-<u'  -L  )t 

I  L  ((!>'  ) 

c 


(5.1) 


Then  we  insert  the  inverse  integral  for  Lc(u')  and  perform  the  u>*  Integration, 
which  yields  the  alternative  expression  for  r(u) 


9 


Ha|k>  *  Jiu)jjk>  ,  (5.3) 

which  defines  the  states  |k>  and  the  eigenvalues  In  turn,  can  be 

represented  with  respect  to  its  own  eigenstates  as 

H.  - 1  •  <5-4) 

k 

where  ■  |k><k|  is  the  projector  onto  |k>.  Then  we  can  express  L^,  Eq.  (3.3), 
in  terms  of  projectors  and  calculate  the  exponential  exp(-iLat).  Finally  we 
obtain 

exp(-iLat)p  »  ^  exp(-iAktT)PkpPt  (5.5) 

kl 

in  terms  of  the  transition  frequencies 


Akt  "  “k  '  wl  ' 

Next  we  substitute  the  definition  (3.5)  of  L£(t)  into  Eq.  (5.2)  and  apply 
the  expansion  (5.5),  which  gives 

r(w)p  -  i  £  Ls(PtSpPkC(Akt-Hj)  -  PkpSPtC(Akt-u>)*)  .  (5.7) 

kt 

This  result  shows  that  r(u>)  can  be  expressed  entirely  in  the  adatom  operator  S 
and  the  Fourier-Laplace  transform  of  the  reservoir  correlation  function  G(t).  In 
order  to  find  the  matrix  representation  of  r(u>),  we  expand  Lg  as  a  commutator 


and  then  insert  the  closure  relation 


in  various  places  in  Eq.  (5.7).  In  terms  of  the  matrix  elements 


Skl  -  <k|S|t>  -  S*k  ,  (5.9) 

we  then  find 

<a|(rU)p)|b>  -  i  £  {sakskteubk+w)<t|p|b> 

kfc 

+  SbkSkiS(Aak'w)*<alpll>} 

-  1  l  (S.kStb5<4U+“> 

kt 

+  S*iSkae(Akb'Ui)*)<k|p|t>  ,  (5.10) 

which  relates  the  matrix  elements  of  T(u>)p  to  the  matrix  elements  of  p.  We 
notice  that  the  frequency  dependence  of  r(w)  only  enters  as  a  shift  of  the 
resonances  A^  in  the  arguments  of  the  correlation  function. 

From  Eq.  (5.4)  we  immediately  find 

<a| ((u-La)p) |b>  «  (w-Aab)<a|p|b>  (5.11) 

for  the  matrix  representation  of  w-L^.  Combination  with  Eq.  (5.10)  then  gives 
the  expansion  of  w-La+ir,  which  is  the  matrix  to  be  inverted  for  evaluation  of 
Po(u>)  from  Eq.  (4.5)  and,  as  we  shall  see  in  due  course,  for  the  evaluation  of 
the  absorption  profile. 

VI.  MARKOV  APPROXIMATION 

In  order  to  illuminate  the  physical  significance  of  a  frequency-dependent 


relaxation  operator,  and  to  establish  the  relation  to  earlier  approaches,  we 


summarize  the  comtnonly-applied  Markov  approximation.  The  phonon  field  has  a 
broad  spectrum  of  closely-spaced  modes  and  a  relatively  high  cut-off  frequency 
which  equals  the  Debye  frequency  for  a  crystal.  These  features  guarantee  that 
the  reservoir  correlation  function  G(t)  decays  to  zero  quite  rapidly  for 
increasing  t,  and  typically  on  a  time  scale  *.  Besides  the  other  two 

typical  time  scales  in  the  problem  are  the  inverse  of  a  level  separation  and 
the  inverse  of  a  relaxation  constant  a^,  which  is  a  matrix  element  of  r(w). 

One  then  asserts  that 

akt  <<  “d  *  (6*° 

for  all  ki,  which  implies  that  Lc(t)  in  the  integrand  of  Eq.  (3.2)  only  deviates 
from  zero  for  t  <  Since  this  is  small  in  comparison  with  the  relaxation 

time,  by  virtue  of  Eq.  (6.1),  we  can  replace  the  time  evolution  of  pQ(t-T)  over 
this  small  time  interval  by  its  free  evolution.  Hence  we  approximate  pQ(t-T)  in 
Eq.  (3.2)  as 

PQ(t-T)  -  exp(iLax)pQ(t)  ,  (6.2) 

and  subsequently  we  can  take  PQ(t)  outside  the  integral.  Then  we  assume  that  we 

are  not  interested  in  the  time  evolution  of  p  (t)  in  a  time  interval  of  the  order 

o 

of  Wp*  after  the  preparation  of  the  initial  state  pQ(0),  which  implies  that  we 
can  replace  the  upper  integration  limit  t  by  infinity.  In  spectral  terms  the 
restriction  (6.1)  means  that  the  width  of  a  line  is  much  smaller  than  the  Debye 
frequency,  and  the  second  approximation  states  that  we  do  not  consider  the 
details  of  the  line  in  the  far  wings,  where  "far"  means  of  the  order  of  from 
the  line  center.  Combining  everything  then  yields  the  equation  of  motion  in  the 
Markov  approximation 


».<*>  ■  <V‘ rM>’0(t> 


(6.3) 


where  the  relaxation  operator  is  now  defined  as 


.  /»  -iL  t  iL 

PM  *  2ii  ^"S  Jn  *  *  Lc(’>*  * 


(6.4) 


The  solution  of  Eq.  (6.3)  reads 


PQ(t)  ■  exp(-i(La*irM)t)po(0)  , 


(6.5) 


which  has  the  Fourier-Laplace  transform 


P0(«)  "  „  i  1T7~  P„(0)  . 

o  uj-L  o 

a  n 


(6.6) 


Comparison  with  our  result  (4.5)  shows  that  the  frequency  dependence  of  the 
relaxation  operator  has  vanished,  and  that  the  Markov  approximation  effectively 
amounts  to  the  substitution  r(w)  -*  ru.  Conversely,  it  then  follows  from  the 
discussion  in  the  first  paragraph  of  this  section  that  the  frequency  dependence 
of  r(w)  reflects  a  time  resolution  on  a  time  scale  of  the  order  of  .  This,  i 
turn,  gives  the  process  a  memory  time  of  the  order  of  w^.  In  the  frequency 
domain  this  has  the  implication  that  now  also  the  line  wings  are  properly 
described,  rather  than  only  the  line  center. 


VI I .  RELATIONS  BETWEEN  Kto)  AND  r. 


From  the  representation  (5.2)  of  r(w)  we  directly  deduce 


|**  dw  exp(-iu)t)r(w)  ■  Lgexp(-iL#T)Lc(t)  ,  t  >  0  , 


(7.1) 


which  reveals  that  the  Fourier-Laplace  inverse  of  T(w)  is  proportional  to  Lc(t), 
and  therefore  it  decays  to  zero  as  a  function  of  t  on  a  time  scale  Then 

T(w)  must  have  a  frequency  width  of  the  order  of  The  right-hand  side  of  Eq. 

(7.1)  also  appears  in  the  definition  (6.4)  of  TM,  which  gives  rise  to  the 


relation 


between  and  r(u>).  Carrying  out  the  time  integral  gives  alternatively 


r 


M 


du  r(u>) 


_i _ 

L  -u  * 
a 


(7.3) 


showing  that  does  not  equal  T(0). 

From  Eq.  (5.11)  we  observe  that  La-w  is  diagonal  with  respect  to  the 
tetradic  adatom  states  |a><b| .  Therefore,  expansion  of  the  inverse  of  L  -w  in 

d 

terms  of  projectors  reads 


_1 _ 

L  -oj 
a 


P 


kl 


kt 


■w  k  1 


(7.4) 


in  terms  of  its  action  on  an  arbitrary  p.  With  the  general  relation  for  Fourier- 
Laplace  transforms 


du'  ;  p((i)'  ) 

UJ-U) 


p(«)  » 


(7.5) 


we  can  evaluate  the  integral  in  Eq.  (7.3),  which  gives 

rMp  *  \  r(Akt)(PkpPl)  *  (7.6) 

kt 

This  simple  relation  between  and  T(w)  displays  that  the  relaxation  operator  in 

the  Markov  approximation  is  determined  by  r(u>),  where  w  only  attains  values  which 

are  equal  to  the  level  separations  A^. 

With  Eq.  (5.10)  we  can  calculate  the  action  of  r(Akt)  on  pkPpt*  We  replace 

P  in  Eq.  (5.10)  by  P^pP^  and  use  the  orthonormality  relation  <a|b>  ■  6^  for  any 

two  states  |a>  and  |b>.  Then  we  sum  over  k  and  t  and  recall  that  A  .  +  A.  "A 

'  '  ab  be  ac 

for  every  a,b,c,  which  gives  the  expansion  of  in  matrix  elements.  We  obtain 


nfimnn 


im an v.'wtwnw wwinn™ ^■■■■■- ^,  JI rj ^  v 


;irj7  mwif  i  .r  i' 
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<aj(rMp)|b>  -  i  £  {Sak  sktC(AtkXt|pjb> 

kl 

+  S*k  SktS(Atk)*<ajpjt» 

•  i  l  <s«kstb«*iJ 

ki 

+  SbtSk»C<4lb)*)<kl|>lt>  •  <7-7> 

12 

which  Is  the  usual  result.  Comparison  with  the  expansion  (5.10)  of  T(u)  shows 
that  the  different  matrix  elements  of  p  are  connected  by  the  same  matrix 
elements  of  S  in  both  cases.  The  only  distinction  is  that  the  correlation 
function  S  appears  with  a  different  argument  in  the  corresponding  matrix  elements 
of  rM  and  r(w).  Since  knowledge  of  the  function  G(u>)  is  already  required  in  the 
Markov  approximation,  no  additional  information  about  the  crystal  or  the  atomic 
states  is  necessary  in  order  to  go  beyond  the  Markov  approximation. 

VIII.  RESERVOIR  CORRELATION  FUNCTION 

Whether  a  Markov  approximation  can  be  justified  or  not  depends  on  the 
behavior  of  the  reservoir  correlation  function  G(t)  from  Eq.  (3.6).  In  order  to 
study  this  issue  quantitatively,  we  adopt  a  Debye  model  for  the  distribution  of 
phonon  modes  in  the  crystal.  Then  the  dispersion  relation  reads 

*s(k)  -  c'kH(wg-c'k)  ,  (8.1) 

with  c*  the  speed  of  sound  and  H  the  unit  step  function.  For  the  thermal- 
equilibrium  density  operator  of  the  crystal  at  temperature  T  we  take 

5  -  (Tr  exp(-H  /kRT))_1  exp(-H  /k_T)  ,  (8.2) 

P  P  P  B  p  B 

where  kg  is  Boltzmann's  constant.  With  the  explicit  form  of  the  reservoir 

o 

operator  R,  Eq.  (2.5),  we  then  immediately  find 


(8.3) 


G(t)  ■  C  dw  u  {(n(w)+l)exp(-iwx)  +  n(w)exp( iwx )}  , 

J(\ 


where  n(w)  equals  the  average  number  of  phonons  in  a  mode  with  frequency  w, 
which  is 


n(w)  ■  {expOiw/kgT)-!} 


(8.4) 


The  parameter  (  is  found  to  be 


x  * 
JiMuip 


(8.5) 


which  depends  only  on  the  mass  of  a  substrate  atom  and  the  cut-off  frequency 
of  the  dispersion  relation. 

It  is  advantageous  to  subdivide  G(x)  into  a  spontaneous  and  a  stimulated 
part  according  to 


G(x)  -  G(t)  +  G(t)  , 
sp  St 


(8.6) 


with  G(t)sjj  by  definition  the  correlation  function  for  T  »  0.  Since  n(w)  «  0  for 
T  *  0,  we  find  with  Eq.  (8.3) 


g(t,  .  fj 

sp  U  /  ^ 

(V) 


c<T).t  ■  «  c 


dw  n(w)ucos(wr )  , 


(8.7) 


(8.8) 


showing  that  the  stimulated  part  is  real,  wheras  the  spontaneous  contribution 
acquires  an  imaginary  component.  In  Figs.  1  and  2  the  correlation  functions 
G(t)gp  and  G(i:)st  are  plotted,  and  it  appears  that  they  indeed  decay  to  zero  on  a 
time  scale  of  the  order  of  w^  However,  the  spontaneous  part  only  disappears 
as  1/(u>jjT),  which  is  very  slow.  This  implies  that  in  integrals  like  in  Eq. 


(3.2),  there  is  a  relatively  large  contribution  from  the  tail  of  G(t),  which 
makes  a  Markov  approximation  at  least  doubtful  when  the  Debye  frequency  becomes 
of  the  order  of  other  frequencies  in  the  problem. 

Matrix  elements  of  the  relaxation  operator  r(u)  are  determined  by  the 
Fourier-Laplace  transform  of  G(t),  which  we  split  up  in  the  same  way.  From  Eq. 
(8.7)  and  Eq.  (8.8)  we  derive 

C(u)sp  *  Curfi(w)H(ujj-u»)  -  iir  ^(wp+tdlogl  l-ojp/u>| )  ,  (8.9) 

2(u»)st  *  C|w|n( |w|  )H(uip-  Jwj ) 

-ix'^PpW  n(u)'  )  — ,  (8.10) 

JO  «'  -w 

where  P  stands  for  principal  value.  The  real  part  of  C(w)  is  only  nonvanishing 

sp 

for  positive  frequencies,  which  has  the  significance  that  the  spontaneous  part  of 
the  relaxation  only  amounts  to  a  decay,  whereas  from  the  relation 

5(-w)st  *  5(w)s^  *  (8.11) 

it  follows  that  the  rate  constants  for  stimulated  excitation  and  decay  of  the 
adatom  are  equal.  In  taking  the  limit  u  ■*  0  in  Eqs.  (8.9)  and  (8.10)  we  obtain 


^  dT  G(T)>p  -  . 

f  dt  G(T)st  -  xCkjT/K  , 


(8.12) 

(8.13) 


which  can  be  regarded  as  a  measure  of  the  relative  contribution  of  the 

spontaneous  and  stimulated  component  to  the  reservoir  correlation  function.  The 

functions  2(<i>)  and  6(w)  are  drawn  in  Figs.  3  and  A,  respectively, 
sp  sc 
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IX.  ABSORPTION  SPECTRUM 

Evaluation  of  the  absorption  spectrum  I(us^)  from  Eq.  (l.S)  is  now 
straightforward.  First  we  introduce  the  correlation  operator 


D(t)  ■  e  iLT[y,p]  , 


(9.1) 


which  is  an  operator  in  the  Hilbert  space  of  the  adatom  and  the  crystal.  Its 
reduced  atomic  part  is  denoted  by  Do(t),  end  the  Fourier-Laplace  transform  CQ(«>) 
is  defined  by  Eq.  (4.1).  Then  Eq.  (l.S)  is  equivalent  to 


I(V  ■  Vt  rfc  Ivsj  Re 


(9.2) 


Next  we  notice  that  D(t)  obeys  the  equation  of  motion 


i  2-  D(t )  -  LD(t )  . 


(9.3) 


which  is  identical  to  Eq.  (2.6)  for  the  complete  density  operator  p(t).  Hence  we 
can  solve  Eq.  (9.3)  for  the  reduced  correlation  operator  along  exactly  the  same 
lines  that  led  to  Eq.  (4.5).  This  gives 


D  (w)  -  .JJrrr}  c  D  (0)  . 

O  Ul-L  +ir(ta>)  o 


(9.4) 


and  the  initial  value  equals 


D  (0)  ■  Tr  lu.p)  -  [p.pj  , 
op  o 


(9.5) 


in  terms  of  the  long-time  solution  pQ  of  Eq.  (4.5).  Combining  everything  then 
yields 


I(uL)  *  ILwL 


fj IT  Re  Tr.  M  VLa+If^)lM*Po1 


(9.6) 


as  a  formal  expression  for  the  absorption  spectrum.  With  the  matrix 
representation  of  T(w)  from  Section  V,  we  can  perform  the  operator  inversion  and 
subsequently  evaluate  I(w^)  for  any  configuration  of  bound  states.  We  notice 
that  the  relaxation  operator  r  in  the  denominator  has  the  incident  frequency 
as  its  argument.  This  reflects  that  the  dependence  of  the  absorption  will 
reveal  the  (non-Markovian)  frequency  dependence  of  the  relaxation  operator. 

Expression  (9.6)  for  I(u)^)  involves  the  steady-state  solution  pQ  of  the 
adatom  density  operator.  From  the  general  identity 

lim  p  (t)  *  p  ■  lim  -iirwp  (w)  ,  (9.7) 

o  o  o 

t-*«®  ur* o 

and  Eq.  (4.5),  it  follows  that  p  is  the  solution  of 

o 

(La-iT(o))po  «  0  .  (9.8) 

Therefore,  it  is  not  necessary  to  calculate  the  Fourier-Laplace  inverse  of  Eq. 
(4.5)  and  then  take  the  limit  t-"»,  which  would  be  a  cumbersome  procedure.  We 
reunark  that  pq  is  completely  determined  by  the  relaxation  operator  at  the  single 
frequency  w  «  0,  in  contrast  to  PQ(t)  for  t  <  <*,  which  involves  T(w)  at  all 
frequencies  u>. 

X.  LINE  SHAPE 

An  adsorbate  potential  V(s)  will  in  general  have  many  bound  states  (-  25), 
and  every  transition  frequency  gives  rise  to  a  spectral  line  around  ■  A (  , 
In  order  to  disentangle  the  contributions  from  the  various  resonances  to  the 
absorption  spectrum  I(u^)>  we  now  consider  the  situation  of  two  levels  |2>  and 
|l>  which  are  separated  by  u.-w.  ■  u  >0.  The  resulting  profile  I(w.  )  is  then 

a!  O  L 


call  led  a  spectral  line. 


Diagonal  matrix  elements  of  the  derivative  of  any  potential  vanish 

identically  if  the  wave  functions  are  eigenfunctions  of  a  Hamiltonian  of  the  form 
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(2.1).  This  implies  ■  S22  *  0.  and  it  follows  from  Eq.  (5.10)  that  the 
relaxation  operator  r(u>)  does  not  couple  between  populations  and  coherences. 

Since  for  two  levels  the  phase  of  a  wave  function  can  have  no  significance,  we 
can  take  the  matrix  elements  and  t0  be  real.**  Then  Eq.  (9.6)  can  be 
rewritten  as 

««L>  ■  Wv52>;  "•  1  <* I ■ L -l  iiroT)*  1 : 2><1 1  - •  I ■ I ■ » I ■ »>  .  C10.1) 

kt 

where  n^  -  <k|pQ|k>  is  the  steady-state  population  of  level  |k>,  and  B  is  the 
Einstein  coefficient  for  stimulated  transition,  defined  as 

B  -  ir(coh2c)"1|M12*cL|2  .  (10.2) 


The  summation  in  Eq.  (10.1)  runs  over  (k,l)  *  (1,2)  and  (2,1)  only. 

Solving  Eq.  (9.8)  for  pq  in  the  case  of  two  levels  is  trivial.  We  only 
mention  that  the  population  difference  Hj-Hj,  which  occurs  in  I(u>^),  becomes 

Re  C(u  ) 

Vn2  *  Re(C(«  )  +2G(«J  )  '  (10,3) 

o  sp  o  St 

We  note  that  this  quantity  is  independent  of  the  interaction  matrix  element  S^ 
and  completely  determined  by  the  reservoir  correlation  function  C(wQ). 
Obviously,  this  is  an  artifact  of  a  two-level  system. 

With  the  basis  |2><l|  and  | 1 ><2 |  we  find 


for  the  relevant  part  of  this  matrix.  Here  the  parameter  function  l(u^)  is 
defined  by 

n(wL)  "  *  S12(C(wL)  +  S(‘WL)*)  *  (10-5) 

which  represents  the  operator  r(uj^)  in  Eq.  (10.4).  Then  the  absorption  line 
shape  is  readily  found  to  be 

2 

4w  w-Re  n(wr  ) 

I(ul  )  *  Ijhif  B(n.-n_)  — = — = - 5 - 5 - j  .  (10.6) 

(wo-uj£+2u)^Imn(u^))  +4u^(Re  n(w^)) 

In  the  Markov  approximation  r)(u».  )  assumes  the  constant  value  n(w  ),  independent 

Li  O 

of  u).  .  Then  n(<*>  )  is  merely  a  parameter,  which  can  be  adjusted  to  fit  a  line 
shape.  It  appears  that  in  our  more  thorough  approach,  the  function  n(“»)  should 
be  evaluated  at  the  frequency  of  the  probe  field.  Especially  the  property 

Re  h(w^)  *  0  ,  ui^  >  Up  ,  (10.7) 


for  a  Debye  model  shows  that  the  absorption  is  identically  zero  for  >  u»p, 
whereas  in  the  Markov  approximation  the  absorption  is  finite  in  this  case.  That 


I(u^)  should  vanish  for  follows  from  energy  conservation. 


a  _  _  1- _ l.  _  J 


photon  must  eventually  end  up  as  an  excitation  of  the  phonon  field,  but  for 


“L  >  “d 


there  are  no  phonon  modes  which  can  accomodate  the  quantum. 


Different 


features  of  the  line  shape  are  illustrated  in  Figs.  5  and  6. 


iw mis  f 1 m  uw  iw  IWBPW-VW  w  wmmtwwwiimmmmmTwwiwiw  w  w»  **  *w*  jt+w*  r* 


XI.  CONCLUSIONS 

Vibrational  relaxation  of  physisorbed  atoms  on  a  harmonic  crystal  is 
studied,  without  the  usual  Markov  approximation.  It  has  been  shown  that  the 
Markov  approximation  can  be  poor  if  the  Debye  frequency  has  an  order  of  magnitude 
which  is  comparable  to  the  widths  and  positions  of  the  spectral  lines,  or  in 
general  where  the  line  wings  are  concerned.  We  have  derived  a  frequency- 
dependent  relaxation  operator  r(u>),  and  from  Eq.  (7.1)  it  follows  that  r(u>) 
deviates  from  a  constant  value  due  to  the  finite  correlation  time  of  the 
reservoir.  This  in  turn  gives  rise  to  a  memory  in  the  time  evolution  of  pQ(t), 
which  is  expressed  by  the  fact  that  pQ(t)  obeys  an  integro-differential  equation 
(Eq.  (3.2))  rather  than  a  first-order  differential  equation.  The  memory  is 
brought  about  by  the  finite  time-width  of  the  operator  Lc(t)  in  the  integrand. 

We  have  solved  the  equation  of  motion  with  a  Fourier-Laplace  transform  and 
applied  the  same  technique  to  evaluate  the  absorption  spectrum.  Significant 
differences  between  the  present  theory  and  its  Markovian  equivalent  are  found  in 
the  spectral  line  profile. 

ACKNOWLEDGMENTS 

This  research  was  supported  by  the  Office  of  Naval  Research  and  the  Air 
Force  Office  of  Scientific  Research  (AFSC),  United  States  Air  Force,  under 
Contract  F49620-86-C-0009.  The  United  States  Government  is  authorized  to 
reproduce  and  distribute  reprints  for  governmental  purposes  notwithstanding  any 
copyright  notation  hereon. 


REFERENCES 


1.  G.  Nienhuis,  Physica  66,  245  (1973). 

2.  S.  Efrima,  K.  F.  Freed,  C.  Jedrzejek  and  H.  Metiu,  Chem.  Phys.  Lett.  74,  43 
(1980). 

3.  Z.  W.  Gortel,  H.  J.  Kreuzer  and  R.  Teshima,  Phys.  Rev.  B  22,  5655  (1980). 

4.  C.  Jedrzejek,  K.  F.  Freed,  S.  Efrima  and  H.  Metiu,  Chem.  Phys.  Lett.  79,  227 
(1981). 

5.  Z.  W.  Gortel,  H.  J.  Kreuzer,  R.  Teshima  and  L.  A.  Turski,  Phys.  Rev.  B  24, 
4456  (1981). 

6.  S.  Efrima,  C.  Jedrzejek,  K.  F.  Freed,  E.  Hood  and  H.  Metiu,  J.  Chem.  Phys. 
79,  2436  (1983). 

7.  A.  I.  Volokitin,  0.  M.  Braun  and  V.  M.  Yakovlev,  Surf.  Sci.  172,  31  (1986). 


8.  A.  A.  Maradudin,  E.  W.  Montroll,  G.  H.  Weiss  and  I.  P.  Ipatova,  Theory  of 
Lattice  Dynamics  in  the  Harmonic  Approximation,  Suppl.  3  of  Solid  State 
Physics  (Academic,  New  York,  1971). 


W.  H.  Louisell, 
York,  1973). 


jantum  Statistical  Properties  of  Radiation  (Wiley,  New 


10.  C.  Cohen-Tannoudji,  in  Frontiers  in  Laser  Spectroscopy,  Proceedings  of  the 
27th  Les  Houches  Summer  School,  edited  by  R.  Balian,  S.  Haroche  and  S. 
Liberman  (North-Holland,  Amsterdam,  1977). 


11.  B.  Fain,  Physica  101A,  67  (1980). 

12.  H.  F.  Arnoldus  and  T.  F.  George,  Phys.  Rev.  B  (1987),  in  press. 

13.  E.  Merzbacher,  Quantum  Mechanics  (Wiley,  New  York,  1961),  p.  41. 

14.  H.  F.  Arnoldus,  S.  van  Smaalen  and  T.  F.  George,  Adv.  Chem.  Phys.  (1987),  in 
press. 


APPENDIX.  DERIVATION  OF  EO.  (3.2)  FROM  E 


In  this  Appendix  we  derive  Eq.  (3.2)  from  the  general  equation  of  motion 
(2.6).  As  an  abbreviation  we  introduce  the  interaction  Liouvillian  as 


LiP  -  -#  1 [ RS , p )  . 


Then  Eq.  (2.6)  reads 


1  j£p(t)  *  <La  +  Lp  +  Lt)p(t)  , 


which  has  to  be  transformed  into  an  appropriate  equation  for  pQ(t).  Integration 


of  Eq.  (A2)  gives 


-i(L  +L  )t 

p(t)  »  e  a  P  pC 0 ) 


ft  *i(L  +L  )(t-t* ).  /.i\ 

■if  df  e  a  P  Lip(t  } 


Next  we  replace  t  by  t'  and  substitute  the  right-hand  side  of  Eq.  (A3)  for  p(t') 
in  the  integrand.  Differentiating  the  result  with  respect  to  time  and  changing 
the  integration  variable  then  yields 


-i(L  +L  )t 


i-j£p(t)  *  (La+L  )p(t)  +  Lie  *  P  p(0) 


ft  -i(L  +L  h 


L.p(t-r) 


as  an  exact  integral  of  Eq.  (A2). 
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As  initial  state  p(0)  we  take  p(0)  ■  Eq.  (3.7),  and  subsequently  we 

factorire  p(t-x)  in  the  integrand,  as  argued  to  be  reasonable  in  section  III. 

Then  we  take  the  trace  over  the  phonon  states  and  insert  definition  (Al)  for  L^. 
The  second  term  on  the  right-hand  side  vanishes  identically  as  follows  from 


-iL  t 

Tr  (R(e  p  pj) 


Tr  Rp 
P  P 


(A5) 


where  the  right-most  equality  relies  on  the  explicity  form  (2.5)  of  R.  The 

appearance  of  two  factors  in  the  last  term  of  Eq.  (A4)  gives  rise  to  four 

terms,  since  L.  is  a  commutator.  Because  R  and  S  commute,  as  do  L  and  L  ,  we 
l  a  p 

can  rearrange  every  term  in  such  a  way  that  all  phonon  operators  are  to  the 
right,  and  all  adatom  factors  are  to  the  left.  As  a  last  step  we  use  R+  *  R,  the 
cyclic  invariance  for  operators  under  a  trace,  and 


-iLt 
e  p 


e  p  e 


(A6) 


for  any  L  and  p.  Comparison  with  the  definition  (3.5)  of  Lc(x)  then  shows  that 
Eq.  (A4)  can  be  cast  in  the  form  (3.2). 


FIGURE  CAPTIONS 


Fig.  1.  Real  (a)  and  imaginary  (b)  part  of  the  spontaneous  component  of  the 

2 

reservoir  correlation  function  as  a  function  of  uipT  and  divided  by  CiUp. 

L  i&j.  Z  •  Stimulated  part  of  the  reservoir  correlation  function,  which  is  real,  as 

2 

a  function  of  Upt  and  divided  by  £ujp.  Curve  (a)  corresponds  to  tfiUp  ■  kgT.  In 
curve  (b)  Oiwp  *  2kgT)  the  plot  is  truncated  at  oipT  *  10. 

Fig.  3.  Real  (a)  and  imaginary  (b)  part  of  G(w)sp  as  a  function  of  w/wp  and 
divided  by  For  u  <  0  and  w  >  Up  the  real  part  vanishes  identically.  The 

singularity  of  the  imaginary  part  at  ui  =  Wp  is  a  result  of  the  sharp  cut-off  Wp 
of  the  dispersion  relation.  For  any  smooth  but  still  arbitrarily  steep  decay  at 


w  *  u»p,  the  value  of  ImG(u)p)sp  would  remain  finite. 

Fig.  4.  Plot  of  the  real  (a)  and  imaginary  (b)  part  of  G(w)st/Cwp  as  a  function 


of  w/wp,  and  for  Itap  *  2kgT.  For  |u|  >  u»p  the  real  part  disappears. 

Line  profile  I(u>p)  as  a  function  of  u)p/wo  anc*  **ivi<*ed  by  ^l^11 
2 

The  parameters  are  Ci»JpS12  *  0.4,  uip  ■  3wq  and  kfiT  «  liuip .  Curve  (b)  is  calculated 


with  Eq.  (10.6)  and  curve  (a)  is  the  Markov  approximation.  The  line  shapes  are 


not  Lorentzians,  as  follows  from  their  asymmetry  around  the  dotted  line  at  uip  * 

w  .  For  (u.  >  u  the  exact  value  and  the  Markov  approximation  always  yield  the 
o  L  o 

sane  value  for  the  absorption. 

■  0.3kgT.  In  this  case  of  relatively  large 
temperature  the  Markov  approximation  (curve  (a))  deviates  considerably  from  the 
exact  value  in  the  blue  wing  (o»p  >  u>o)  of  the  line.  Especially  for  (Up  >  ujp  the 
Markov  approximation  predicts  a  finite  absorption,  which  cannot  be  correct. 


Same  as  Fig.  5  but  with 
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